Limnology and Oceanography

The Canadian Arctic Archipelago (CAA) is warming at an unprecedented rate, leading to sea ice loss and glacial retreat. Marine-terminating (tidewater) glaciers can fuel summertime marine productivity by delivering nutrient-rich deep waters via upwelling to the surface ocean. While the impact of glacier-induced upwelling has been well-studied in the context of phytoplankton and primary productivity, its effects on broader marine microbial communities remain poorly understood. We investigated how glacier-driven upwelling shapes marine microbial (bacterial and archaeal) communities across a series of sites in the CAA. At upwelling sites, the upper 50 m of the water column exhibited elevated nutrient concentrations and physical characteristics that resembled deeper waters, which were associated with differences in microbial community composition relative to non-upwelling sites. Our results indicate that upwelling influences microbial communities in surface waters in two ways. It directly introduces typically deeper-water-associated taxa into surface waters and reshapes ecological niches by enhancing nutrient supply and stimulating primary production, indirectly driving changes in microbial communities. The enrichment of Candidatus Nitrosopumilus, a deep water nitrifier, likely affects nitrogen cycling and raises the possibility of active nitrification in surface waters. Likewise, the increased abundance of taxa known to be associated with phytoplankton-derived organic matter in upwelling regions suggests an enhanced capacity to process organic matter generated from elevated primary productivity. Ultimately, as tidewater glaciers continue to retreat, the resulting changes in the glacially-driven upwelling regime will likely shift marine microbial communities towards assemblages adapted to less productive ecosystems, with implications for nutrient cycling in these systems. ImportanceClimate change has a disproportionate impact on the Arctic, with rising temperatures causing increased marine-terminating glacier retreat and changes in the marine water column structure. The consequent loss of the ability of these glaciers to upwell deep water to the surface ocean results in a reduction of nutrient delivery and mixing in these ecosystems. Previous work has highlighted the importance of marine-terminating glaciers in sustaining phytoplankton productivity during the summer season through this delivery of deep-water nutrients to the surface ocean. The impact of glacially-induced upwelling on marine bacterial and archaeal communities, however, remains underexplored. We found that in regions with glacially-driven upwelling, the surface ocean showed enrichment of phytoplankton-associated taxa and nitrifiers commonly associated with deep waters. This work underscores the role of glacially-driven upwelling in structuring both microbial communities and nutrient cycling, suggesting that glacier loss could reshape community composition and biogeochemical processes in a rapidly changing Arctic.

Ocean warming is altering abiotic environments and biotic interactions experienced by marine organisms, where sensitive early developmental windows occur in biologically complex seawater communities. The impact of these interactions on developmental processes and fitness in hosts is not well understood, but likely contingent on the establishment of a host-associated microbiome. Here, we hypothesize that temperature and microbial exposure during embryogenesis influence larval microbiome assembly and host morphology. Strongylocentrotus purpuratus embryos were raised in low microbial richness (LMR) or high microbial richness (HMR) seawater at ambient (14 {degrees}C) or elevated (18 {degrees}C) temperature, then collected at 2, 4, and 6 days post-fertilization (dpf) following multiple feedings. Higher microbial diversity was observed in larvae that developed in HMR seawater when compared to LMR. Differences in relative abundances of dominant microbial families between seawater and larvae suggest some degree of host selectivity in microbiome assembly. Temperature did not strongly alter microbiome composition, but both temperature and microbial condition led to differences in larval morphology by 6 dpf, potentially due to enrichment of microbes with chemoheterotrophic functions. By linking how temperature and microbial communities interact with host development, we contribute novel insights into how early-life environmental conditions impact holobiont formation and morphology. One sentence summaryEarly developmental temperature and microbial conditions shape larval microbiome establishment and morphology.

Near-shore coral reefs in southern Guam (Mariana Islands) experience severe sedimentation, in particular during the wet season when rainfall and erosion are high. We sampled fragments of the reef-forming coral Porites lobata from opposite ends of a sedimentation gradient in Fouha Bay, southern Guam, during dry and wet seasons. Using DNA metabarcoding, we characterized the diversity and composition of P. lobata-associated Symbiodiniaceae and bacterial microbiome communities. As in many species of Porites, Symbiodiniaceae communities of P. lobata were dominated by variants of Cladocopium C15 with sites showing differences in Symbiodiniaceae communities attributable to variation in these Cladocopium C15 variants. Bacterial microbiomes of P. lobata were dominated by Endozoicomonadaceae, a family of putative coral bacterial endosymbionts involved in nutrient cycling. Site and seasonal differences in bacterial diversity and community composition were apparent. In close proximity to the mouth of the river draining into Fouha Bay, bacterial diversity was highest during the wet season when sedimentation is generally severe. Microbiome reorganization in response to sedimentation may explain this result, but we also found overrepresentation of bacteria associated with terrestrial origin close to the river mouth and/or during the wet season. Together these patterns highlight that coral Symbiodiniaceae and bacterial communities are both spatially and temporally structured in this disturbed system. IMPORTANCEThis study provides a time series dataset of coral-associated microorganisms, including dinoflagellate algae and bacteria, from a tropical bay impacted by sedimentation that results from upstream erosion of disturbed soils. Characterizing temporal patterns of coral-associated microbes provides insights into the dynamic nature of these communities. While microbiome variability across sites and seasons may be a result of acclimatization to different environmental conditions, we identified bacterial groups of putative terrestrial origin in sampled coral microbiomes that may have been exported from eroded soils to the near-shore reef. Considering that disturbed soils act as hotspots for the proliferation of potentially harmful substances, such as antimicrobial resistance genes, understanding microbial community connections at the marine-freshwater-terrestrial interface is an important step toward evaluating environmental impacts across connected ecosystems from ridge to reef.

Environmental DNA (eDNA) surveys are increasingly used to assess marine biodiversity and inform deep-sea environmental decision-making, including mineral resource management and fisheries oversight. Yet standard low-volume protocols inherited from coastal work may be inadequate at depth, and no quantitative framework links depth and ecosystem context to defensible filtration volume targets. We compiled 841 eDNA samples from eight expeditions across the North Atlantic, Wider Caribbean, and Pacific (surface to 4000 m) to quantify how recoverable eDNA scales with depth and surface productivity, and to derive depth- and productivity-aware sampling targets. Total eDNA concentration declined with depth as a power law, with attenuation exponents (b) modulated by surface productivity: most gradual in eutrophic waters (b = 0.67), intermediate in mesotrophic (b = 0.90), and steepest in oligotrophic systems (b = 1.25); volume-weighted models explained 66-88% of the variance. At a fixed extract-concentration target, required filtration volumes diverged ~7-fold between oligotrophic and eutrophic systems at 200 m and ~38-fold at 4000 m. Conventional Niskin sampling, therefore, undersamples deep-sea biodiversity, particularly in mid- to low-productivity systems. Among laboratory parameters, the assay-specific extract-concentration target exerted greater leverage on required volume than extraction efficiency or elution volume. Volume-aware sampling paired with optimized recovery should be routine in deep-sea eDNA surveys.

The mechanisms driving host-symbiont associations across space and time in contemporary mutualisms can give insight into the capacity for symbiotic organisms to respond to environmental change. High specificity between partners can increase cooperation and facilitate efficient holobiont selection, whereas low specificity could reduce host benefit, but facilitate adaptive associations across heterogeneous environments. The present study explores specificity in natural populations of a cnidarian-algal model, Exaiptasia diaphana, across a latitudinal gradient to understand the genetic and environmental effects driving host-symbiont associations, and their relation to heritable and/or environmental symbiont acquisition. We found that symbiotic associations were extremely flexible in E. diaphana, regardless of transmission mode. E. diaphana were capable of associating with diverse symbiont communities across genetically identical hosts seeded with vertically transmitted symbionts, as well as across highly connected host populations which acquire symbionts horizontally. Host population connectivity was complex and unrelated to geographic distance, whereas symbiont community composition tracked the thermal gradient, potentially due to context dependent biotic interactions. These results indicate that in a flexible symbiosis, symbiont communities are environmentally-determined, suggesting the future of this symbiosis will likely depend on climate adaptation of symbionts.

The prevalence of nitrogen limitation and nitrogen-phosphorus co-limitation (henceforth referred to as nitrogen-related limitation) in Norwegian lakes and their relationships with atmospheric nitrogen deposition, climate, dissolved organic matter (DOM), and catchment characteristics were assessed across space and time. Routine monitoring data from 1,529 lakes in the national Vannmiljo database were analyzed for two multi-year periods (1995-2009 and 2010-2025). Limitation was inferred using the molar NO--N/TP ratio as an indicator of dissolved inorganic nitrogen availability. Nitrogen-related limitation was widespread in both periods and exhibited strong regional structure, with highest prevalence in northern regions and lowest prevalence in southwestern Norway. Overall prevalence increased from 31% to 38% between periods, with significant increases in western regions. Regional-scale models identified climate, forest cover, DOM, agriculture, and atmospheric nitrogen deposition as predictors of limitation probability, whereas study period per se and bog/peatland cover were not significant. At the local scale, atmospheric nitrogen deposition and DOM were the only consistent predictors, with substantially lower explanatory power than at the regional scale. These results indicate that large-scale environmental gradients play a major role in shaping nutrient stoichiometry in Norwegian lakes. Because the monitoring dataset primarily represents lakes affected by human activities, the findings are particularly relevant for water management. The widespread occurrence of nitrogen-related limitation suggests that nitrogen availability may influence phytoplankton growth in many systems and that dual-nutrient management strategies addressing both nitrogen and phosphorus may be required in many regions.

Quantifying the lipid biosynthesis rate of archaea in hot spring sediments is necessary to interpret the abundance, isotopic patterns, and environmental significance of archaeal lipid biosignatures, with implications for modern biogeochemical cycling and astrobiology. Here, we performed lipid hydrogen stable isotope probing (LH-SIP) experiments on whole sediments collected from two high-temperature, suboxic, circumneutral hot springs in Yellowstone National Park (USA) and El Tatio Geyserfield (Chile). We determined the incorporation of 2H2O into intact polar lipids (IPLs) which provides a taxon- and metabolism-agnostic quantification of biosynthesis under near-natural conditions. We targeted isoprenoid glycerol dialkyl glycerol tetraether lipids (IPL iGDGTs) and recovered structures with 0 to 7 cyclopentyl rings from both springs. We observed minor 2H-uptake into archaeal IPLs in spring sediments in Yellowstone, corresponding to decadal-scale apparent generation times (16 {+/-} 7 years), and no uptake in El Tatio sediments (consistent with minimum generation times of 35 {+/-} 5 years). We infer that net production of sedimentary IPL-iGDGTs is very slow, consistent with a combination of slow archaeal growth, persistence of older IPLs, lipid recycling, and/or contributions from recently sedimented planktonic biomass. These are the first direct, ex situ estimates of archaeal lipid production rates in terrestrial hydrothermal systems using LH-SIP incubations and provide critical constraints for interpreting archaeal lipids in ancient hot spring deposits. This research establishes a framework for assessing activity by slow-growing extremophilic archaea in hydrothermal environments and provides support for targeting hydrothermal deposits on Mars for biosignature detection efforts. Plain Language SummaryHot springs on Earth are important natural laboratories for understanding how signs of life might form and be preserved in hydrothermal environments on early Earth or Mars. In this study, we examine the rate of archaeal lipid biosignature production in sediments from two hot springs in Yellowstone National Park and the El Tatio Geyserfield in Chile. We used a method that measures new microbial production by tracing heavy hydrogen from labeled water as microbes incorporate that hydrogen into newly made lipids in their cell membranes. We found that archaeal lipids in hot spring sediments are produced very slowly, on timescales of decades. This result, along with the chemical stability of lipids and the rapid mineralization rate in hot springs, may allow these molecular biosignatures to be entombed and preserved in hot spring mineral deposits. These results help us better interpret ancient hydrothermal deposits on Earth and support the idea that slowly growing microbial communities could still leave detectable molecular traces in similar environments on Mars and other rocky planets. Key PointsO_LILipid hydrogen stable isotope probing is applied to high temperature hot spring sediments for the first time C_LIO_LIIn hot spring sediments, archaeal lipid production occurs on decadal timescales comparable to some marine sediments, but are much faster than the century- to millennia-scale rates observed in the deep subsurface C_LIO_LIConfirmation of archaeal lipid synthesis in hot spring sediments adds additional support for targeting Martian hydrothermal deposits for biosignature detection efforts C_LI

O_LIThe frequency and intensity of heat events is increasing across marine and terrestrial ecosystems. Within the same ecological community, the relative exposure and sensitivity to heat stress may vary considerably among interacting species, like predators and prey. This can be especially true for species that interact at the aquatic-terrestrial interface, as well as for interactions between primarily nocturnal and diurnal species, making it difficult to predict how such communities will respond to habitat warming. C_LIO_LIThermal limit metrics such as CTmax are often assumed to equate with ecological death because such temperatures impair behavioral activity and/or physiological functioning. Prey that are diurnally active can be more frequently exposed to temperatures that approach CTmax compared to their nocturnal predators, which may use thermal refuges during the day. Yet the impacts of daytime heat exposure on nighttime predation risk remain unknown. C_LIO_LIHere, we compared the thermal environment, performance, and heat tolerance between the predatory blue crab, Callinectus sapidus and one of its prey species, the mangrove periwinkle Littoraria anguilifera in a tropical mangrove ecosystem. We examined how exposing prey to heat stress at and below their CTmax affected their capacity to avoid predation in the field at night when predation risk is highest. C_LIO_LIWe found that acute exposure to temperatures near CTmax during the day increased the prey species susceptibility to predation during recovery at night. Although both interacting predator and prey have high thermal tolerance, prey are exposed to conditions that already reach CTmax, suggesting that current extremes in temperatures may already be influencing vulnerability to predation in this ecosystem. C_LIO_LIOur results suggest that differential exposure to sublethal heat stress in diurnal prey relative to their predator, along with the subsequent impact of these exposures on predation risk, will play a role in shaping these interacting as climate warms. C_LI

Reef restoration practitioners aim to preserve coral genetic diversity by protecting reefs and cultivating diverse genotypes in coral nurseries. However, cryptic genetic lineages in most corals complicate restoration strategies, as the role of between-lineage genetic divergence remains unclear regarding adaptation. In Montastraea cavernosa, researchers have identified cryptic lineages, some strongly segregated by depth. We conducted a ten-week reciprocal transplantation experiment using two cryptic lineages restricted to shallow water (<10m depth), with one lineage more common on nearshore reefs and the other on offshore reefs. We aimed to quantify lineage-specific responses to the environment that explain the genetic and ecological divergence between the two lineages. Surprisingly, the strongest response was not lineage-specific. Instead, both lineages exhibited strong and similar changes in growth and metabolomic profiles, depending on the transplantation habitat. These results suggest that cryptic lineages employ similar mechanisms of adaptation and acclimatization to environmental challenges, despite their genetic distinction.

Linking animal movements to environmental drivers is essential for understanding ecological processes and anticipating species responses to climate change. We investigated habitat-specific movements in a globally significant aggregation of loggerhead turtles (Caretta caretta) nesting in Cabo Verde. Satellite tags on 15 adults (12 females, 3 males) provided multi-year tracks spanning breeding, migration, and foraging habitats. Movements and phenology differed by habitat. During the breeding season, females used either coastal areas, remaining within [~]20 m depth, or undertook long looping forays up to 360 km. Males showed two strategies: two remained resident in Cabo Verde waters, including Fra, the largest male tracked (Curved carapace length of 105 cm compared with a male mean of 90.7 {+/-} 10.3 cm), while the third migrated annually to distant foraging grounds and returned ahead of the subsequent breeding season. In foraging habitats, turtles adopted neritic or oceanic strategies: neritic turtles remained localised in warm, productive waters, whereas oceanic turtles ranged widely in deeper, less productive areas. Time- and space-shift analyses showed that oceanic foragers used intermediate sea surface temperature and chlorophyll-a conditions relative to nearby or temporally shifted alternatives, consistent with movement within a thermal-trophic trade-off. Together, these results show how sex, body size, and energy balance drive habitat-specific movement dynamics in a changing ocean.

Microplastics (MP) are widespread in aquatic ecosystems and pose a threat to freshwater biodiversity. While numerous studies examine physiological effects on aquatic organisms, less is known about how MP alter chemically mediated interactions that regulate predator-prey dynamics. Predator-induced defenses in Daphnia depend on detecting kairomones and represent an important form of adaptive phenotypic plasticity. Whether MP interfere with these responses, and through which mechanisms, remains unclear. Here, we show that polystyrene MP impair predator-induced defenses across Daphnia species by disrupting predator-cue-mediated plasticity at the behavioral, morphological, and molecular levels. In D. longicephala, chronic exposure to PS fragments weakened Notonecta-induced morphological defenses, whereas additive-containing PS fragments nearly suppressed defense formation and reduced body size. Consistent with these phenotypic effects, proteomic analyses revealed alterations in pathways related to molting and chitin metabolism, linking MP exposure to impaired defense formation. In D. magna, PS particles attenuated fish kairomone-induced diel vertical migration, with stronger effects for larger particles, consistent with reduced effective availability or perception of predator cues. Natural limestone particles caused only minor effects, indicating particle-specific rather than general particle-driven responses. Our findings demonstrate that MP can disrupt adaptive predator-prey interactions with potential cascading consequences for freshwater food webs.

Guinea worm disease eradication efforts may benefit from environmental surveillance methods capable of detecting infected copepod intermediate hosts in aquatic habitats. We developed a three-dimensional, spatially explicit agent-based model to examine how ecological processes influence detection probability for a hypothetical water sampling method. The results show that surveillance sensitivity is shaped by the combined effects of larval diffusion, copepod density, and pond size, with interactions among these factors producing nonlinear relationships. Detection, in our model, was concentrated within a relatively restricted period after larvae matured to the infective stage and before dispersal and mortality reduced presence, indicating a limited spatiotemporal window for effective sampling. Surveillance performance peaked under intermediate dispersal regimes that generated sufficient spatial overlap between larvae and intermediate hosts, while both limited dispersal and excessive diffusion reduced detection by constraining encounters or diluting larval concentrations. Increasing habitat size reduced detection by diluting larval concentrations, but the magnitude of this effect depended on copepod density and dispersal dynamics, producing nonlinear and threshold responses rather than simple scaling with pond volume. Spatial and temporal patterns of detection shifted as larvae dispersed, with the most favorable detection periods occurring when both larval abundance and intermediate host encounters were elevated. These findings indicate that surveillance can be guided by local ecological conditions. When the timing of larval introduction is uncertain, effective surveillance requires repeated sampling over time to capture transient windows of detectability and the sampling will be less effective in very stagnant and highly mixed waterbodies. Overall, this study demonstrates how mechanistic modeling can support the design and interpretation of environmental surveillance strategies for Guinea worm eradication programs. Author summaryGuinea worm disease is close to eradication but confirming that transmission has fully stopped remains difficult because detecting infectious larvae in water is challenging. Transmission depends on freshwater copepods that become infected after ingesting Guinea worm larvae. These copepods are short-lived and unevenly distributed within ponds, and infected individuals may die before larvae reach the infective stage. As a result, environmental detection is inherently uncertain. We developed a three-dimensional agent-based model to simulate larval dispersal, copepod infection, and water sampling in a pond environment. The model shows that detection is constrained to a brief period when mature larvae and copepods overlap in space and time, and that this window depends strongly on local ecological conditions such as larval dispersal, copepod density, and pond size. Because infected copepods can be present outside these narrow detection windows, negative water samples do not necessarily indicate absence of transmission, highlighting the need for repeated, spatially targeted surveillance during the final stages of eradication.

Coral reefs deliver vital services via a complex three-dimensional framework sustained by the balance between calcium carbonate production and erosion, or the net carbonate budget state. In many tropical western Atlantic reefs, ecological decline has reduced carbonate production, yielding near-neutral or negative budgets. Yet some reefs retain high coral cover and, theoretically, should also have high net positive budgets, yet often show modest carbonate accumulation. We used the remote reef of Cayo Arenas in the Campeche Bank, Gulf of Mexico, to test whether in reefs under suboptimal (variable) environmental conditions, high coral production is offset by robust bioeroder communities, producing neutral budgets. At 14 sites, we quantified carbonate producers and bioeroders to estimate gross production, bioerosion, and net budget states. Despite relatively high live coral cover, mean net carbonate budgets were approximately neutral. Crucially, this neutrality arose not from depressed biological activity (as in degraded reefs) but from an active equilibrium: vigorous carbonate production coupled with substantial bioerosion. These reefs, therefore, represent a contemporary, functional reef state in net stasis. Distinguishing active-neutral from impoverishment-neutral regimes is critical for predicting reef trajectories under environmental change and for targeting management, although near-stasis emerging from high carbonate turnover can appear functionally intact yet operate with limited buffering capacity against net carbonate loss.

Global changes in land use and nutrient cycling are transforming ecosystems at unprecedented rates, with significant consequences for infectious disease dynamics. Aquatic environments are particularly vulnerable because the interplay of habitat modification, nutrient enrichment, and biodiversity loss can drive pronounced changes in the community composition of food webs, including hosts and parasites. Yet, despite well-documented effects of habitat modification on aquatic communities and food webs, the mechanisms through which these changes influence infectious disease dynamics remain poorly resolved. This gap arises, in part, because it remains challenging to disentangle how multiple stressors interact to shape disease outcomes and quantify parasite levels and host densities from field-collected samples. Here, we illustrate two tools that might help address these challenges. First, highly sensitive droplet digital PCR can quantify infection loads even when the signal:noise ratio is low. Second, stepwise Bayesian path analyses can identify the direct and indirect pathways connecting land-use changes to infectious disease dynamics. As a case study, we examined cyclopoid copepods and their helminth parasite, Schistocephalus solidus, across 47 freshwater lakes on Vancouver Island, a region strongly shaped by commercial logging, including widespread clear-cutting of old-growth forests. Our results reveal a positive correlation between copepod density and deforestation, potentially mediated by associated changes in water quality and calanoid copepods, key competitors of the focal host. ddPCR enabled sensitive detection of extremely low parasite signals in field-collected copepods. We detected positive infections in only 19.5% of the lakes surveyed, highlighting the difficulty of assessing disease dynamics in natural populations. Nonetheless, this study highlights the challenges of linking land-use change to disease outcomes, while also demonstrating that sensitive molecular and statistical tools offer new ways to reveal these hidden connections.

Many environmental bacteria do not readily grow under laboratory conditions and population establishment often occurs stochastically. Although the scout hypothesis has been proposed to explain stochastic population establishment in environmental bacteria, how stochastic population establishment is shaped by individual cell growth behaviors in environmental isolates remains unclear. In the present study, we focused on the ammonia-oxidizing bacterium Nitrosomonas sp. PY1 and showed that environmentally responsive individual cell growth behavior, incorporating time-dependent stochastic growth initiation, shapes both deterministic and stochastic population establishment dynamics. Using single-cell observation, we revealed that PY1 altered cell growth behavior in response to surrounding biomass production ({Delta}Vt). These {Delta}Vt-dependent changes in growth behavior were suppressed by the addition of its own cell-free supernatant (CFS), indicating the presence of a growth regulation mechanism via cell-cell communication. Replicate cultures under the same conditions showed that the population establishment of PY1 was stochastic, whereas the model strain Nitrosomonas europaea exhibited synchronized population establishment, consistent with previous reports. This stochasticity in PY1 was also eliminated by the addition of CFS. Finally, a simulation model based on {Delta}Vt-dependent cell growth behavior of PY1 successfully reproduced synchronized population establishment in the presence of CFS. By contrast, the stochastic population establishment observed in the absence of CFS was successfully reproduced by a model incorporating {Delta}Vt-independent growth initiation following a Weibull distribution. Such environmentally responsive changes in population establishment dynamics may contribute to the low isolation success of environmental bacteria and sudden blooms of the rare biosphere.

The maintenance of endosymbiosis in cnidarians depends on the tight regulation of host immunity, cell cycle, and nutrient exchange, yet how these processes are impacted by interacting environmental stressors remains largely unknown. To address this, we employed physiological metrics, gene expression analysis, microbiome characterization, imaging (NF-{kappa}B localization, endoplasmic reticulum ultrastructure, EdU labeling), and stable isotope tracing in the model sea anemone Exaiptasia diaphana to examine the effects of heat and nitrate on these regulatory processes, individually and in combination. Heat treatment led to NF-{kappa}B activation, proteostatic stress, suppression of nutrient exchange, decreased cell-cycle progression, and microbiome restructuring, with all effects more pronounced in symbiotic than aposymbiotic anemones. In symbiotic anemones, nitrate partially offset these heat-induced responses through sustained carbon translocation, suggesting that the presence of symbionts, in conjunction with elevated nitrate, can temporarily buffer host thermal stress. However, prolonged combined exposure resulted in holobiont failure. These findings reveal that while nitrate enrichment can transiently delay the onset of bleaching, it does not preserve the regulatory networks required for symbiotic stability -- underscoring the vulnerability of cnidarian holobionts to the compounding effects of warming and nitrate pollution.

Antimicrobial resistance (AMR) genes are increasingly recognized as an emerging environmental contaminant. Yet, the ecological mechanisms shaping their distribution across natural landscapes remain poorly understood. Here, we quantified AMR gene abundances in microbial communities sampled from wild fish from eight freshwater lakes on Vancouver Island and paired these gene-level measurements with fine-scale limnological and land-use data. Using droplet digital PCR, field surveys, and an iterative spatial forecasting framework that integrates Random Forest models with regression kriging, we explored how watershed-scale processes relate to variation in AMR genes across lakes. Our analyses reveal potential associations between elevated AMR gene levels, changes in water quality, deforestation, and geographic proximity to salmon aquaculture. By integrating data across biological and spatial scales, from genes within microbial communities to lake-level conditions and landscape patterns, this study illustrates the value of combining quantitative molecular measurements with geospatial modeling to identify environmental factors that may promote antimicrobial resistance in natural systems. Our approach provides a proof-of-concept and a general predictive framework for generating hypotheses and informing future monitoring efforts aimed at understanding, managing, and forecasting environmental reservoirs of resistance. SignificanceAntimicrobial resistance (AMR) genes are ancient components of environmental microbiomes. Yet, the mechanisms that generate modern hotspots of resistance across natural landscapes remain unclear. Here, we reveal how watershed-scale environmental change, including water quality metrics linked with deforestation and proximity to salmon aquaculture, predicts elevated AMR gene levels in the microbiomes of wild fish populations. By combining quantitative droplet digital PCR with ecological data and geospatial modeling, we move beyond isolated surveillance data to identify ecological mechanisms that promote antimicrobial resistance in freshwater ecosystems. This integrative approach provides mechanistic insight into why certain habitats, and the organisms within them, become reservoirs of resistance while others do not. Our findings highlight the importance of ecological context in understanding resistance evolution and offer a predictive tool for informing proactive monitoring and management strategies.

Plankton are essential to marine ecosystems, supporting food webs and mediating biogeochemical processes such as carbon export to depth. Their spatial distribution influences ecosystem dynamics and serves as an indicator of environmental change. Although drifting plankton could be expected to exhibit random distribution, numerous studies have revealed significant heterogeneity in their spatial patterns. However, very few studies targeted plankton distribution at the centimeter scale in situ, despite its importance for understanding biological processes. We argue that centimeter-scale distances in plankton could reveal potential ecological interactions. Using an extensive in situ dataset of 18 million planktonic organisms collected by the In Situ Ichthyoplankton Imaging System (ISIIS), which images multiple organisms simultaneously and preserves their positions in the water column, we analyzed centimeter-scale distances in plankton. By comparing observed distances with those expected under a random distribution, we assessed potential interactions at three levels: among all organisms, within plankton groups and across groups. Our results show that planktonic organisms exhibit non-random distributions at the centimeter scale, with smaller distances than expected, suggesting potential ecological interactions. Notably, distances up to 11 cm were the most informative, which is much larger than typical interaction distances in plankton. Additionally, observed distances were compatible with a simple attraction model. Finally, we propose the non-randomness of distances as a novel metric of interaction strength in plankton ecological networks and compare it against classical empirical or co-occurrence networks. These results offer new insights into in situ interactions and how they shape plankton distribution at centimeter scale. Significance statementThis study reveals that planktonic organisms exhibit non-random spatial distributions at the centimeter scale, highlighting the importance of ecological interactions in shaping their distribution at this scale. By analyzing an extensive in situ plankton imaging dataset, we introduce a novel metric of interaction strength based on the non-randomness of distances between organisms, and compare it to common interaction metrics. These findings challenge the traditional view of plankton as passive drifters by highlighting that their distribution at microscale is shaped not just by physical processes such as turbulence but also by ecological interactions. Author contributionsJOI contributed to data acquisition. TP processed the data under the supervision of JOI. TP, JOI and BBC designed the study. TP conducted the analyses under the supervision of MF, JOI and BBC. TP wrote the initial draft of the manuscript. All authors contributed to the interpretation of results, supported manuscript preparation and approved the final submitted version.

The fish skin microbiome serves as a protective barrier, influencing host health and facilitating interactions between the host and its environment. While several studies have characterised the composition and roles of the fish skin microbiome, there remains a paucity of data on how environmental variation influences these microbes in natural populations. Here, we used 16S rRNA gene sequencing to characterise the skin microbiome of wild three-spined stickleback populations and examine how environmental factors influence microbial diversity and community composition across 17 freshwater lochs on the island of North Uist, Scotland. Analysis of 239 samples revealed a set of dominant bacterial genera commonly associated with other fish species, including Janthinobacterium, Pseudomonas, Acinetobacter, and Psychrobacter, that constituted a core skin microbiota across lochs. Microbiome composition was primarily shaped by environmental variables, particularly habitat, water pH, conductivity, and metal concentrations, with pH emerging as a key driver of community structure. Host sex also influenced microbiome variation, with several taxa differing in relative abundance between males and females. Alpha-diversity was higher among stickleback fish from lochs with a neutral pH compared with those from alkaline and acidic environments. Differential abundance analyses identified 27 and 24 amplicon sequence variants (ASVs), respectfully, associated with variations in pH and host sex, including members of Psychrobacter, Sphingobacterium, Carnobacterium, Chryseobacterium, and Arthrobacter, highlighting the combined influence of environmental and host factors on microbiome composition in wild fish populations in freshwater environments.

Understanding risks of biological invasions associated with ballast water (BW) requires full understanding of the biodiversity transported in ballast tanks. Here we characterize the remarkable level of diversity that can be carried in the BW of a single vessel. To maximize our ability to capture BW diversity we: 1) utilized DNA-based methods to describe biodiversity, including both native and non-native taxa; 2) exploited multiple primer sets targeting multiple genomic loci with different expectations for taxonomic coverage; 3) assessed multiple tanks on a single vessel to capture different communities present in different tanks; and 4) sampled those tanks with far higher-than-usual replication both to improve representation of the diversity present and to enable statistical estimation of total richness. Using this approach, we found extraordinarily high diversity associated with a single vessel. Across all loci, we estimate a total of 272 taxa that can be assigned species names; looking more broadly at unnamed molecular operational taxonomic units, our estimates are between 425 and 742 individual taxa, depending on the locus. We confirm that only a fraction of this diversity would be captured with typical sampling efforts. We found that different loci capture different snapshots of biodiversity and that our ability to detect taxa of interest (e.g., non-native species) depends on sampling effort and genomic locus. Our results expand upon previous studies describing highly diverse BW communities and add to a growing literature that demonstrates the value of molecular methods for characterizing those communities and assessing the associated risk of non-native species introduction.